This is an example document for single cell data analysis using RamEx, replicating the results of the relevant analysis in RamEx: An R package for high-throughput microbial ramanome analyses with accurate quality assessment. bioRxiv

RamEx could be downloaded from [here] (https://github.com/qibebt-bioinfo/RamEx)

Test dataset could be download from ScienceDB

Load RamEx

library(RamEx)
library(dplyr)

Data Loading

plot(data)
 [1] "0h"   "2h"   "4h"   "8h"   "12h"  "20h"  "28h"  "36h"  "44h"  "52h"  "60h"  "68h"  "76h"  "84h"  "92h"  "100h"

Reductions

data.reduction <- Feature.Reduction.Pca(data[, data@wavenumber < 1700 | data@wavenumber > 2700], n_pc = 20) %>% Feature.Reduction.Umap(., min.dist = 0.5, n_neighbors = 50)

Phenotype analysis

clusters <- Phenotype.Analysis.Louvaincluster(data, resolutions = c(0.4), n_pc=20)
Creating the communities...
Processing resolution: 0.4 
Merging the isolate samples...
Plot.reductions(data.reduction, color=clusters$Resolution_0.4)

mean.spec(data$normalized.data, clusters$Resolution_0.4 , gap=0.3)
[1] "1" "2" "3" "4"

Plot.Distribution(data$group, clusters$Resolution_0.4 )

Raman markers

Plot.Heatmap.Markers(data, markers$markers_singular$wave)
Plot.Markers.Spectrum(object = data, markers=markers)

markers_population_cluster <- Raman.Markers.Correlations(data, group = clusters$Resolution_0.4, by_average = T, min.cor = 0.95, extract_num = FALSE,min.range = 20)
Finding singular markers ... 
Plot.Markers.Spectrum(data, markers_population_cluster)

MCR-ALS

mcr <- Spectral.Decomposition.Mcrals(data,n_comp = 2)
Plot.ViolinBox(mcr$concentration, data@meta.data$group)

Plot.ViolinBox(mcr$concentration, clusters$Resolution_0.4)

mean.spec(t(mcr$components), as.factor(colnames(mcr$components)), gap=0)
[1] "Comp 1" "Comp 2"

IRCA

IRCA_interests <- Intraramanome.Analysis.Irca.Global(data)
It may take a while to drawing the plot, please be patient.
It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

It may take a while to drawing the plot, please be patient.

band_ann <- data.frame(Wave_num=c(781,847,886,915,995,1075,1117,1148,1235,1295,1328,1438,1573,1591,1647,1738,2840,2900,2919),Group=c('Nucleic Acid','Polysaccharide','Polysaccharide','Polysaccharide','Protein','Lipids+Protein','Polysaccharide','Protein','Protein','Lipids','Protein','Lipids','Nucleic Acid','Lipids','Lipids','Lipids','Lipids','Lipids','Lipids'))
Intraramanome.Analysis.Irca.Local(data, bands_ann = band_ann)

LS0tDQp0aXRsZTogIlllYXN0cyINCm91dHB1dDogaHRtbF9ub3RlYm9vaw0KLS0tDQoNClRoaXMgaXMgYW4gZXhhbXBsZSBkb2N1bWVudCBmb3Igc2luZ2xlIGNlbGwgZGF0YSBhbmFseXNpcyB1c2luZyBSYW1FeCwgcmVwbGljYXRpbmcgdGhlIHJlc3VsdHMgb2YgdGhlIHJlbGV2YW50IGFuYWx5c2lzIGluIFtSYW1FeDogQW4gUiBwYWNrYWdlIGZvciBoaWdoLXRocm91Z2hwdXQgbWljcm9iaWFsIHJhbWFub21lIGFuYWx5c2VzIHdpdGggYWNjdXJhdGUgcXVhbGl0eSBhc3Nlc3NtZW50XShodHRwczovL2RvaS5vcmcvMTAuMTEwMS8yMDI1LjAzLjEwLjY0MjUwNSkuICpiaW9SeGl2KiANCg0KUmFtRXggY291bGQgYmUgZG93bmxvYWRlZCBmcm9tIFtoZXJlXSAoaHR0cHM6Ly9naXRodWIuY29tL3FpYmVidC1iaW9pbmZvL1JhbUV4KQ0KDQpUZXN0IGRhdGFzZXQgY291bGQgYmUgZG93bmxvYWQgZnJvbSBbU2NpZW5jZURCXShodHRwczovL3d3dy5zY2lkYi5jbi9lbi9maWxlP2ZpZD01ODI3MGJmYzhlNzQxYTM1NmEzNTQ0ZmIzMDEyZDllNiZtb2RlPXN1cmwmc3VybD1WTmZqdTYpDQoNCiMjIExvYWQgUmFtRXgNCmBgYHtyfQ0KbGlicmFyeShSYW1FeCkNCmxpYnJhcnkoZHBseXIpDQpgYGANCg0KIyMgRGF0YSBMb2FkaW5nDQpgYGB7cn0NCmRhdGEgPC0gcmVhZFJEUygiVGVzdF9SYW1FeF9kYXRhX3llYXN0LnJkcyIpDQpwbG90KGRhdGEpDQpgYGANCg0KIyMgUmVkdWN0aW9ucw0KYGBge3J9DQpkYXRhLnJlZHVjdGlvbiA8LSBGZWF0dXJlLlJlZHVjdGlvbi5QY2EoZGF0YVssIGRhdGFAd2F2ZW51bWJlciA8IDE3MDAgfCBkYXRhQHdhdmVudW1iZXIgPiAyNzAwXSwgbl9wYyA9IDIwKSAlPiUgRmVhdHVyZS5SZWR1Y3Rpb24uVW1hcCguLCBtaW4uZGlzdCA9IDAuNSwgbl9uZWlnaGJvcnMgPSA1MCkNCmBgYA0KIyMgUGhlbm90eXBlIGFuYWx5c2lzDQpgYGB7cn0NCmNsdXN0ZXJzIDwtIFBoZW5vdHlwZS5BbmFseXNpcy5Mb3V2YWluY2x1c3RlcihkYXRhLCByZXNvbHV0aW9ucyA9IGMoMC40KSwgbl9wYz0yMCkNClBsb3QucmVkdWN0aW9ucyhkYXRhLnJlZHVjdGlvbiwgY29sb3I9Y2x1c3RlcnMkUmVzb2x1dGlvbl8wLjQpDQptZWFuLnNwZWMoZGF0YSRub3JtYWxpemVkLmRhdGEsIGNsdXN0ZXJzJFJlc29sdXRpb25fMC40ICwgZ2FwPTAuMykNClBsb3QuRGlzdHJpYnV0aW9uKGRhdGEkZ3JvdXAsIGNsdXN0ZXJzJFJlc29sdXRpb25fMC40ICkNCmBgYA0KDQojIyBSYW1hbiBtYXJrZXJzDQpgYGB7cn0NCm1hcmtlcnMgPC0gUmFtYW4uTWFya2Vycy5Db3JyZWxhdGlvbnMoZGF0YS5yZWR1Y3Rpb24sIG1pbi5jb3I9MC43LCBwYWlyZWQgPSBUKQ0KUGxvdC5IZWF0bWFwLk1hcmtlcnMoZGF0YSwgbWFya2VycyRtYXJrZXJzX3Npbmd1bGFyJHdhdmUpDQpQbG90Lk1hcmtlcnMuU3BlY3RydW0ob2JqZWN0ID0gZGF0YSwgbWFya2Vycz1tYXJrZXJzKQ0KbWFya2Vyc19wb3B1bGF0aW9uX2NsdXN0ZXIgPC0gUmFtYW4uTWFya2Vycy5Db3JyZWxhdGlvbnMoZGF0YSwgZ3JvdXAgPSBjbHVzdGVycyRSZXNvbHV0aW9uXzAuNCwgYnlfYXZlcmFnZSA9IFQsIG1pbi5jb3IgPSAwLjk1LCBleHRyYWN0X251bSA9IEZBTFNFLG1pbi5yYW5nZSA9IDIwKQ0KUGxvdC5NYXJrZXJzLlNwZWN0cnVtKGRhdGEsIG1hcmtlcnNfcG9wdWxhdGlvbl9jbHVzdGVyKQ0KYGBgDQoNCiMjIE1DUi1BTFMNCmBgYHtyfQ0KbWNyIDwtIFNwZWN0cmFsLkRlY29tcG9zaXRpb24uTWNyYWxzKGRhdGEsbl9jb21wID0gMikNClBsb3QuVmlvbGluQm94KG1jciRjb25jZW50cmF0aW9uLCBkYXRhQG1ldGEuZGF0YSRncm91cCkNClBsb3QuVmlvbGluQm94KG1jciRjb25jZW50cmF0aW9uLCBjbHVzdGVycyRSZXNvbHV0aW9uXzAuNCkNCm1lYW4uc3BlYyh0KG1jciRjb21wb25lbnRzKSwgYXMuZmFjdG9yKGNvbG5hbWVzKG1jciRjb21wb25lbnRzKSksIGdhcD0wKQ0KYGBgDQoNCiMjIElSQ0ENCmBgYHtyfQ0KSVJDQV9pbnRlcmVzdHMgPC0gSW50cmFyYW1hbm9tZS5BbmFseXNpcy5JcmNhLkdsb2JhbChkYXRhKQ0KYmFuZF9hbm4gPC0gZGF0YS5mcmFtZShXYXZlX251bT1jKDc4MSw4NDcsODg2LDkxNSw5OTUsMTA3NSwxMTE3LDExNDgsMTIzNSwxMjk1LDEzMjgsMTQzOCwxNTczLDE1OTEsMTY0NywxNzM4LDI4NDAsMjkwMCwyOTE5KSxHcm91cD1jKCdOdWNsZWljIEFjaWQnLCdQb2x5c2FjY2hhcmlkZScsJ1BvbHlzYWNjaGFyaWRlJywnUG9seXNhY2NoYXJpZGUnLCdQcm90ZWluJywnTGlwaWRzK1Byb3RlaW4nLCdQb2x5c2FjY2hhcmlkZScsJ1Byb3RlaW4nLCdQcm90ZWluJywnTGlwaWRzJywnUHJvdGVpbicsJ0xpcGlkcycsJ051Y2xlaWMgQWNpZCcsJ0xpcGlkcycsJ0xpcGlkcycsJ0xpcGlkcycsJ0xpcGlkcycsJ0xpcGlkcycsJ0xpcGlkcycpKQ0KSW50cmFyYW1hbm9tZS5BbmFseXNpcy5JcmNhLkxvY2FsKGRhdGEsIGJhbmRzX2FubiA9IGJhbmRfYW5uKQ0KYGBgDQo=